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SUMMARY 



'The  author'-e  method  to  compute  the  steady  low-speed  inviscid  flow  past  a 
wing  in  free  air,  is  extended  to  take  account  of  ground  effect.  The  basic  method 
represents  the  perturbation  due  to  the  wing  by  iteratively  computed  distributions 
of  sources  and  doublets  on  the  wing  chordal  surface;  at  each  iteration  the  ground 
effect  is  represented  by  the  images  in  the  ground  plane  of  these  distributions, 
the  strengths  of  which  are  calculated  from  the  computed  errors  in  the  boundary 
conditions  on  upper  and  lower  surfaces.  For  a given  angle  of  incidence  (and 
Mach  number),  several  heights  above  the  ground  in  succession  can  be  treated,  with 
some  economy  in  computing  time. 


Results  are  presented  for  a two-dimensional  section  (RAE  100)  and  for  two 

• 1. 

variants  of  a wing  of  airbus  type.  Comparisons  with  another  method  for  the 
RAE  100  section  suggest  that  the  present  method  generally  needs  at  least  three 
iterations,  and  that  for  very  accurate  results  at  low  ground  heights  further  work 


on  the  program  is  needed. 
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1 INTRODUCTION 

In  earlier  work*  we  have  shown  how  to  compute  the  steady  inviscid  flow 
past  a thick,  cambered  and  twisted  wing  (neglecting  dihedral)  in  free  air,  at 
low  Mach  number.  The  flow  field  is  represented  by  suitable  source  and  doublet 
(singularity)  distributions  in  the  wing  chordal  surface,  the  strengths  of  which 

are  calculated  iteratively  by  the  second-order  small-perturbation  theory  of 

2 . 3 4.. 

Weber  . A computer  subroutine,  developed  by  Ledger  and  Sells  , is  available 

to  calculate  the  corresponding  flow  fields. 

In  this  approach  to  the  problem  of  a wing  in  free  air,  the  iterative 
calculation  of  the  strengths  of  the  singularity  distributions  is  based  at  each 
step  on  finding  the  amount  by  which  the  velocity  fields  violate  the  wing 
boundary  conditions  of  zero  normal  flow  on  upper  and  lower  surfaces.  For  the 
problem  of  a wing  near  the  ground,  this  idea  can  be  extended  very  simply.  The 
effect  of  a solid  planar  ground  can  be  represented  by  image  singularity  distri- 
butions, and  the  corresponding  interference  velocity  fields  induced  on  the  wing 
can  also  be  calculated  by  the  Ledger-Sells  subroutine  and  added  to  the  basic 
velocity  fields.  The  necessary  modifications  are  set  out  in  section  2.  The 
interference  fields  are  expected  to  be  small  con^iared  with  the  basic  velocity 
fields,  and  the  major  consequence  of  this  is  that  they  should  produce  only  small 
perturbations  in  the  boundary-condition  errors,  which  can  easily  be  absorbed 
into  the  general  error  fields  to  be  cancelled  (approximately)  at  each  iteration. 
Thus,  no  special  modification  is  needed  in  the  program  to  calculate  the 
successive  singularity  strengths. 

Most  of  the  computer  time  in  this  method  is  needed  by  the  Ledger-Sells 
integration  subroutine,  and  at  first  sight  it  might  appear  that  the  additional 
computation  of  the  interference  velocity  fields  on  upper  and  lower  wing  surfaces 
would  increase  run  times  by  a factor  of  about  three.  In  fact  the  penalty  is  not 
as  large  as  this,  for  two  reasons.  First,  since  the  numerical  variation  in  the 
interference-velocity  integrands  subsides  rapidly  with  distance  from  the 
singularity  plane,  the  integration  (which  is  adaptive,  i.e.  compares  results  for 
successively  larger  numbers  of  function  evaluations)  takes  less  time  than  that 
for  the  basic  velocity  fields.  Secondly,  the  differences  between  the  inter- 
ference fields  on  upper  and  lower  surfaces  is  also  expected  to  be  small  compared 
with  the  basic  fields,  and  indeed  may  be  small  compared  with  the  interference 
fields  themselves,  so  that  in  most  circumstances  it  should  suffice  (and  this 
has  been  assumed  in  writing  the  program)  to  calculate  these  interference  fields 


I J 
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on  one  surface  only.  For  very  accurate  work  this  assumption  may  fail,  as  we 
shall  see  in  section  4 (Results),  and  then  it  may  be  necessary  to  determine  the 
separate  upper-  and  lower-surface  interference  fields,  perhaps  using  Taylor 
series  as  in  Ref. I (where  it  was  shown  how  to  tackle  the  similar  problem  for 
cambered  wings) . 

The  method  assumes  the  planar  wake  of  linear  theory,  and  is  therefore 

8 

less  general  than  that  of  Maskew  , who  is  able  to  distribute  his  singularities 
(vortex  rings)  over  the  wing  camber  surface  and  one  of  two  possible  wake  models. 
However,  the  results  reported  by  Maskew  do  suggest  that  the  difference  in 
results,  due  to  different  wake  models,  is  often  small  compared  with  the 
difference  due  to  ground  effect  itself,  so  we  might  argue  that  if  we  could  deter- 
mine the  wake  shape  as  part  of  the  solution,  the  results  would  be  unlikely  to 
differ  much  from  those  obtained  with  a sin^jle  wake  model.  More  significant  may 
be  the  fact  that  we  can  take  the  thickness  interaction  into  account  as  well  as 
those  of  camber,  twist  and  incidence,  and  in  the  results  section  we  show  a 
case  for  which  the  extra  lift  generated  by  ground  effect  is  due  principally  to 
.hanges  in  the  lower-surface  pressure  distribution,  for  which  the  boundary 
layer  is  relatively  thin,  while  the  upper-surface  pressure  distribution  ai  ' its 
associated  thick  boundary  layer  do  not  change  much  from  their  free-air  disposi- 
tions. This  would  imply  that  for  this  case  viscous  effects  are  similar  to  those 
for  the  wing  in  free  air;  a result  which  would  be  unlikely  to  emerge  from  an 
application  of  Maskew's  thin-wing  method.  This  is  an  example  of  the  general 
contention,  which  we  now  make,  that  the  inclusion  of  thickness  interactions 
might  allow  better  prediction  of  viscous  effects. 

A possible  further  advantage  of  the  present  method  is  that,  whereas  (for 
one  configuration)  a panel  and  matrix  method  might  need  up  to  four  times  the 
effort  to  calculate  influence  coefficients  in  the  presence  of  ground  effect, 
compared  with  the  corresponding  free-air  problem,  the  present  method  needs  rather 
less  than  three  times  the  total  effort,  './hen  results  are  required  for  more  than 
one  ground  height,  the  present  method  can  carry  over  the  stored  basic  velocity 
field  from  a previous  iteration,  thus  saving  one  application  of  the  Ledger-Sells 
subroutine;  but  a panel  and  matrix  method  can  also  save  time  by  holding  the  sub- 
matrix  of  self-induced  influence  coefficients  in  store,  and  it  is  not  easy  to 
predict  which  method  would  gain  on  this  point. 

The  program  is  written  in  Fortran  and,  with  arrays  dimensioned  for  a 
12  X 15  collocation  grid,  requires  about  57K  words  of  core  store  (including 
about  13K  words  for  a vortex  lattice  matrix). 
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2 THE  IMAGE  SYSTEM  IN  GROUND  EFFECT 

We  consider  a finite  wing  moving  at  incidence  a and  unit  speed,  parallel 
to  an  infinite  plane  representing  the  ground.  (When  the  wing  is  relatively 
near  the  ground,  departures  from  parallel  motion  must  be  small.)  The  flow  is 
considered  inviscid  and  irrotational , but  may  be  slightly  compressible,  even  at 
landing  and  take-off  speeds.  We  superpose  unit  velocity  on  the  system  in  the 
usual  way,  so  that  the  wing  is  considered  at  rest  in  a moving  airstream,  and 
define  cartesian  coordinates  (x,y,z)  with  origin  T^  at  the  wing  root  trailing- 
edge,  with  the  x-axis  downstream  and  parallel  to  the  ground,  the  y-axis  to  star- 
board and  the  z-axis  vertically  upwards. 

In  Ref. I we  considered  a wing  in  free  air,  with  no  ground  effect,  and 

2 

showed  how  to  calculate  the  resulting  flow  field  iteratively,  following  Weber  . 
The  perturbation  flow  field  due  to  the  wing  is  represented  by  distributions  of 
sources  and  doublets  on  the  wing  chordal  surface.  To  start  the  iteration,  the 
strengths  of  these  singularities  are  estimated  by  linear  theory;  the  resulting 
velocity  fields  are  computed  at  certain  collocation  points  on  the  upper  and  lower 
wing  surfaces,  the  errors  in  the  boundary  conditions  of  zero  normal  flow  are 
determined  and  combined  into  symmetrical  and  antisymmetrical  parts,  R^,  R^^ 
respectively,  and  perturbation  singularity  distributions  are  generated  to  cancel 
these  errors  approximately,  again  using  linear  theory.  The  cycle  can  be 
repeated  as  often  as  desired;  in  practice  two  or  three  iterations  suffice  for 
acceptable  convergence. 

When  the  wing  chordal  surface  is  non-planar,  we  consider  a plane  parallel 
to  the  y-axis,  containing  the  local  chordline  at  the  collocation  wing  section 
y “ constant  and  inclined  to  the  onset  flow  and  the  x-axis  at  an  angle 
a + ».p(y)  « where  a^(y)  is  the  local  twist  angle,  and  assume  that  a singularity 
distribution  on  this  plane  produces  the  same  velocity  field  - to  second-order 
accuracy  - at  this  collocation  section,  as  would  the  same  singularity  distribu- 
tion on  the  twisted  chordal  surface.  Such  an  assumption  has  to  be  made  here, 

since  we  intend  to  compute  the  perturbation  velocity  fields  using  the  Ledger- 
3 4 

Sells  * computer  subroutine  which  applies  only  to  singularity  distributions 
in  a plane. 

In  extending  the  method,  the  effect  of  ground  proximity  can  now  be 
represented  by  planar  image  source  and  doublet  singularities  (of  the  same  sign 
for  sources,  of  opposite  sign  for  doublets),  and  these  can  easily  be  incorporated 
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into  the  iterative  scheme.  As  the  Ledger-Sells  subroutine  can  be  used  at  any 
points  off  the  plane  containing  the  singularities,  we  can  calculate  the  inter- 
ference velocity  fields  induced  on  the  wing  by  the  image  singularities,  as  well 
as  by  the  basic  chordal-surface  singularities,  and  take  them  into  account  when 
we  compute  the  residual  errors  R^,  in  the  boundary  conditions.  We  expect 

that  on  the  wing  surface  the  interference  fields  will  be  small  compared  with 
the  basic  velocity  fields,  so  that  the  residual  errors  will  not  be  grossly 
increased  (i.e.  by  an  order  of  magnitude  or  more)  and  we  can  again  cancel  them 
approximately  by  adding  suitable  perturbation  singularities  as  before,  with 
unimpaired  convergence.  Thus  no  special  modification  is  needed,  even  to  the 
first  set  of  singularity  distributions,  estimated  from  linear  theory. 

A further  consequence  of  the  comparative  separation  of  the  wing  and  its 
image  is  that  the  difference  between  the  interference  velocities  on  the  upper 
and  lower  surfaces  of  a thin  wing  near  a given  point  P on  the  wing  chordal 
surface  will  be  comparatively  small,  and  so  for  moderate  accuracy  it  will 
suffice  to  approximate  to  these  on  both  surfaces  by  the  values  at  one  point 
only,  say  P , with  a saving  in  computer  time.  The  possible  level  of  error  due 
to  this  assumption  is  discussed,  for  a typical  case,  in  section  4. 

To  compute  the  interference  velocities,  we  need  to  know  the  coordinates  of 
the  points  P relative  to  a local  cartesian  coordinate  system  (x',y,z')  with 
origin  in  Che  image  singularity  plane  corresponding  to  the  section  y = constant  . 
We  take  the  x'-axis  parallel  to  the  image  of  the  local  section  chordline,  and  the 
z'-axis  normal  to  it  and  in  the  general  upward  direction.  We  denote  the  ground 
clearance  height  of  the  local  section  trailing-edge  T by  g^(y,a)  , so  that  the 
distance  between  T and  its  image  T'  is  2g^  . An  expression  for  g^(y,a) 
is  derived  in  the  Appendix.  The  leading-edge  image  L'  has  the  coordinates 
(x'  = x^(y) , y,  z'  =0),  and  so  T'  has  coordinates  (x^(y)  + c(y) , y,  0)  . 

If  the  point  P (usually  referred  to  as  a *Weber  point'*)  has  the  local 
percentage-chord  coordinate  C , with  5=0  representing  the  leading-edge  L , 
then  the  distance  PT  is  (I  - 5)c(y)  . Referring  to  Fig.l,  we  can  now  write 
down  the  coordinates  of  P as: 

[xj^  + c + 2g^  sin  a*  - ( 1 - 5)c  cos  2a*,  y,  2g^  cos  a*  + (I  - 5)c  sin  2a*] 
where  a*  = a + a^(y)  . 
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The  total  interference  velocity  u , due  to  sources  and  doublets,  is 

“(j 

computed  as  (u'.v'.w')  in  the  (x',y,z')  system,  and  we  need  the  components 
G b G 

(u„,v  ,w  ) in  a local  (x,y,z)  system  defined  in  the  singularity  plane  containing 
G G G 

the  real  local  chordline,  as  in  Ref.l.  These  are 

u = u'  cos  2a*  - w ' sin  2a* 

G G G 


w_  = u'  sin  2a*  + w'  cos  2a* 
G G G 


Since  the  interference  velocity  components  do  not  change  sign  as  we  cross 
the  wing  chordal  surface,  the  complete  velocity  field  now  has  components  (U,V,W) 
in  the  local  (x,y,z)  system,  given  on  upper  and  lower  surfaces  (upper  and  lower 
signs  respectively)  by 


u - * .(,)  ± Q, 

V - (Qj  * v^)  = Qj 
W = (Q^  + w^)  ± 

where  ^ are  defined  in  equations  (3-7)  to  (3-12)  of  Ref.l  for  the  wing  out 
of  ground  effect. 

3 SUCCESSIVE  CALCULATION  FOR  SEVERAL  AIRCRAFT  HEIGHTS 

It  is  convenient  to  be  able  to  make  calculations  for  several  aircraft 
heights  in  one  computer  run.  If  the  incidence  and  Mach  number  are  held  constant, 
the  perturbation  source  and  doublet  distributions  should  not  vary  greatly  as  the 
height  is  varied,  and  so  it  should  be  practicable  to  use  the  distributions  found 
iteratively  for  one  height  as  a starting  point  for  the  next  value  of  height.  The 
corresponding  velocity  fields  generated  directly  on  the  wing  are  still  in  store, 
and  so  the  only  preparatory  work  needed  is  the  recalculation , at  the  new  relative 
positions  of  wing  and  ground  image,  of  the  interference  velocity  field  u^  due 
to  the  total  source  and  doublet  distributions  corresponding  to  the  previous 
height,  which  must  of  course  be  stored. 
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In  this  way  we  can  save  the  major  part  of  one  iteration.  So,  having 
chosen  (as  input  data)  a number  of  iterations  judged  necessary  to  secure 

acceptable  accuracy,  it  seems  reasonable  to  reduce  by  one  after  completing 

the  calculations  for  the  first  value  of  height,  and  this  is  now  done  by  the 
program. 

If  the  ratio  of  height  above  ground  to  semi -span  is  sufficiently  large, 
the  problem  reduces  to  that  of  an  isolated  wing,  and  no  calculations  of  u 
need  be  made.  In  the  program,  the  value  100  of  this  ratio  is  judged 
sufficiently  large. 

4 RESULTS 

4 . 1 RAE  100  section 

As  a first  check  on  the  method,  we  compare  lift  coefficients  in  two- 
dimensional  flow,  calculated  for  the  10%  thick  RAE  100  section,  with  some 
results  obtained^  by  the  panel  method  of  Hess  and  Smith  (AMOS)^.  The  section 

K 1 

is  given  by 


~ = 0.148l88v'C(l  - O (l  - I 0 . ^<0.75, 


Z 

— = 0.0855564(1  - O , C > 0.75. 

c 


To  obtain  a nearly  two-dimensional  flow  without  rewriting  the  program,  this 

section  was  embodied  in  a rectangular  wing  with  aspect  ratio  60.  Twelve  chord- 

wise  collocation  points  were  taken.  For  the  two  incidences  a = ±5°  the 

section  lift  coefficient  is  plotted  against  the  reciprocal  of  the  ratio  of 

the  height  above  ground  h^  (of  the  quarter  chord  point)  to  section  chord  c , 

in  Fig. 3.  With  this  choice  of  abscissa  we  can  include  the  results  for  the 

isolated  wing,  c/h^  = 0 , for  which  case  the  two  methods  agree  well.  For 

finite  h , differences  appear;  as  h /c  decreases,  the  difference  in  C due 

W W L 

to  ground  effect  is  at  first  overpredicted,  then  later  underpredicted,  by  the 
present  method  relative  to  the  AMOS  method,  the  cross-over  occurring  at  roughly 
hp/c  * 0.4  . 

There  are  three  possible  reasons  for  these  discrepancies.  The  first  is 
that  the  difference  in  due  to  ground  effect  may  be  relatively  sensitive  to 
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the  level  of  accuracy  with  which  the  boundary  conditions  are  satisfied.  For  the 
cases  h^/c  = 0.75  and  0.25,  both  at  5°  incidence,  two  sets  of  results  are 
available,  obtained  using  different  first  guesses  while  testing  the  program's 
option  to  deal  with  several  heights  consecutively.  In  the  table  below  we  show 
the  numerically  largest  residual  errors  R^,  in  the  boundary  conditions,  and 

the  corresponding  values  obtained  for  . (The  results  with  the  smaller 

residual  errors  are  the  ones  shown  in  Fig. 3.)  The  number  of  iterations  for 
these  runs  was  = 2 . 


u 

O' 

42 

max 

max 

'^L 

0.75 

-0.0009 

0.0000 

AM 

0.0018 

0.0003 

OS 

0.6122 

0.6155 

0.6007 

0.25 

-0.0068 

-0.0042 

AM 

0.0092 

0.0054 

3S 

0.6537 

0.6793 

0.7058 

As  R and  R decrease,  at  h /c  = 0.75  , C converges  away  from  the 

u X>  L* 

AMOS  value,  but  at  =0.25  it  converges  towards  the  AMOS  value.  In  this 

latter  case,  extrapolating  roughly  to  R^  = 0 gives  = 0.7206  while  extra- 
polating to  R = 0 gives  C = 0.7156  , so  it  is  reasonable  to  assume  that 

X/  Li 

the  final  converged  value  lies  somewhere  in  this  range.  In  any  event,  the  size 
of  these  differences  for  this  level  of  accuracy  in  the  boundary  conditions, 
which  would  be  quite  reasonable  for  the  wing  out  of  ground  effect,  suggests  that 
it  may  be  worthwhile  to  ask  for  at  least  = 3 iterations  in  general. 

Another  possible  source  of  error  is  the  approximation  of  equal  interference 
velocities  on  upper  and  lower  wing  surfaces;  to  see  by  how  much  this  approxima- 
tion might  be  in  error,  the  author  computed,  from  two-dimensional  formulae,  the 
interference  velocity  fields  due  to  the  basic  source  and  doublet  distributions: 


= 2 


dz 
t 

cdx 


= 4a 


(1  - o 
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He  found  that  for  the  smallest  height  ratio  = 0.25  and  5°  incidence, 

the  difference  in  interference  streamwash  u on  the  upper  and  lower  surfaces 
was  typically  about  25%  of  the  mean  interference  streamwash  and  about  2%  of 
the  total  streamwash.  Uhich  of  these  two  numbers  is  more  relevant  is  difficult 
to  say,  but  (looking  again  at  the  table  and  comparing  with  the  value  0.5810  of 
for  the  section  out  of  ground  effect)  we  see  that  the  error  could  be  about 
10%  of  the  difference  due  to  interference,  and  about  1%  of  the  total.  So  it 
may  be  that  for  very  accurate  work  the  labour-saving  assumption  of  equal  inter- 
ference velocities  on  upper  and  lower  surfaces  will  have  to  be  discarded,  and 
the  program  revised  accordingly. 

The  third  possible  source  of  error  is  just  the  discretization  error  due 
to  the  finite  number  of  collocation  points  in  either  solution.  However,  the 
good  agreement  shown  in  Fig. 3 for  the  section  out  of  ground  effect  suggests 
that  this  error  may  be  comparatively  small,  though  it  is  unwise  to  draw  any 
definite  conclusion  thereon  without  comparing  the  actual  pressure  distributions. 

4 . 2 Finite  wings 

To  demonstrate  the  program's  performance  on  a finite  wing,  we  present  some 
results  for  two  variants  of  a wing  of  'airbus'  type,  the  planform  of  which  is 
shown  in  Fig. 4.  To  get  some  idea  of  what  happens  at  a wing  tip  comparatively 
near  the  ground  one  variant  (denoted  Var.l)  is  so  chosen  that,  at  incidence 
-0.05°,  the  trailing-edge  is  everywhere  at  constant  height  g.^,  above  the  ground. 
When  this  constant  height  is  about  half  the  wing  span  (g.p/s  = 1.0),  there  is 
little  difference  from  the  isolated  wing  condition.  Figs. 4 to  6 show  what  happens 
when  g.j,/s  is  decreased  to  0.25  and  then  to  0.10.  In  this  range,  at  the  root 
and  in  mid-semispan  the  section  lift  (Fig. 4)  is  increasing  and  we  observe  from 
the  chordwise  plots  (Figs.5a-c)  that  the  increase  is  due  mainly  to  a change  in 
the  lower  surface  pressure  distribution,  the  upper  surface  distribucion  changing 
comparatively  little,  so  that  (for  this  case)  viscous  effects  due  to  boundary 
layer  changes  are  likely  to  be  small.  This  could  not  easily  be  predicted  by  a 
method  for  a thin  wing,  which  takes  no  account  of  the  interaction  of  thickness 

g 

and  ground  effect,  such  as  Maskew's  *.  Near  the  tip,  however,  the  section  lift 
starts  to  decrease  somewhere  between  g.j,/s  = 0.25  and  0.10;  there  is  clearly 
a further  three-dimensional  effect  outboard,  and  the  detailed  computer  output 
shows  that  the  induced  incidence  (due  to  interference)  is  decreasing  in  this 

* In  the  present  method,  even  though  the  perturbation  due  to  the  ground  is  taken 
to  be  the  same  on  both  surfaces,  the  effect  on  Cp  of  compounding  the 
velocity  vectors  is  not  the  same. 


k 


I 


1 1 


region,  as  the  ground  clearance  is  reduced.  The  local  spanwise  centre  of 
pressure  plot  (Fig. 6)  tells  us  also  that  the  overall  centre  of  pressure  moves 
forward  as  the  ground  clearance  decreases,  which  could  imply  an  unfavourable 
nose-up  pitching  moment  (for  this  hypothetical  case). 

To  demonstrate  the  ability  of  the  program  to  handle  wings  in  which  the 
trai 1 ing-edge  ground  clearance  is  not  constant  along  the  span,  we  consider 
another  variant  (denoted  Var.2)  of  this  wing  in  which  the  ground  clearance  ratio 
(at  a = -0.05°)  varies  linearly  from  0.10  at  the  root  to  0.05  at  the  tip.  In 
Figs. 7 and  8 the  spanwise  lift  distribution  and  centre. of  pressure  locus  are 
shown  for  Var.2,  and  also  for  Var.l  with  g.^/s  = 0.10  and  with  g.^/s  = 0.05  , 
as  before.  As  one  might  expect,  the  results  for  Var.2  agree  with  the  results 
for  Var.l  at  g^/ s = 0.10  inboard,  and  with  the  other  set  outboard.  Incident- 
ally, the  trend  towards  loss  of  lift  with  decreasing  ground  clearance  has  spread 

inboard  to  the  root,  as  we  expect  to  happen  in  due  course  from  two-dimensional 
9 

studies  (Kiichemann  ) . 

No  quantitative  conclusions  are  drawn  from  these  results,  for  two  reasons. 
Firstly,  the  problem  is  somewhat  academic;  secondly,  the  differences  shown  might 
be  comparable  with  the  computational  errors  discussed  earlier  with  the  results 
for  the  two-dimensional  section,  though  if  these  errors  are  assumed  to  be  similar 
for  all  runs,  which  in  this  case  all  started  from  the  same  initial  conditions, 
the  qualitative  comparisons  above  are  still  likely  to  be  sensible. 

5 CO.NCLUSION 

The  method  of  chordal-surface  singularities  for  calculating  flow  past  an 

2 1 

isolated  wing,  proposed  by  Weber  and  implemented  by  Sells  , has  been  extended 
to  cover  the  case  of  a wing  in  ground  effect. 

By  taking  advantage  of  the  comparatively  small  size  of  the  interference 
velocity  fields,  we  have  avoided  extensive  changes  to  the  basic  program;  more- 
over, computing  times  are  not  greatly  increased,  though  results  for  a two- 
dimensional  section  suggest  that  the  pressure  distribution  and  section  lift  may 
be  more  sensitive  to  tiie  residual  errors  in  the  wing  boundary  conditions  than 
in  the  isolated-wing  problem,  so  that  at  least  one  more  iteration  may  be 
advisable.  Since  the  basic  program  requires  at  least  two  iterations,  this 
implies  that  for  the  same  accuracy  the  new  program  ought  to  be  run  for  at  least 
three  . 


The  inclusion  of  the  effect  of  wing  thickness,  as  well  as  those  of  camber, 
twist  and  incidence,  on  the  surface  pressure  distributions,  by  allowing  a 
better  prediction  of  boundary-layer  effects,  may  serve  to  fill  a gap  in  current 
numerical  technique  for  this  problem.  As  an  example  of  this  hypothesis,  we  show 
results  for  a case  in  which  we  can  deduce  that  viscous  effects  are  probably 
similar  to  those  for  the  wing  in  free  air. 


Appendix 

THE  LOCAL  TRAIL DiG-EDGE  GROUND  CLEARANCE 
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In  order  to  compute  the  interference  velocities  induced  on  the  wing  by 
the  image  singularities,  we  need  some  information  about  the  ground  clearance. 
Since  the  wing  trailing-edge  is  frequently  made  up  of  straight  line  segments, 
and  is  generally  nearer  than  the  leading-edge  to  the  ground,  it  seems  reasonable 
to  work  with  trailing-edge  ordinates. 

We  input,  as  part  of  the  flow  field  data,  the  height  (ground  clearance) 
h.j  of  the  section  trailing-edge  Tq  at  the  root  y - 0 . As  part  of  the  data 
specifying  the  wing  geometry  we  input  the  spanwise  variation  z^(y)  of  the 
height  of  the  local  trailing-edge  T^  above  the  root  trailing-edge  T^  in  the 
reference  configuration  of  zero  incidence,  o ■ 0 ; see  Fig. 2,  in  which  T^  has 
been  projected  onto  the  root  plane  y “ 0 , which  is  the  plane  of  the  paper. 
Knowing  the  local  leading-edge  ordinate  *^(y)  ®nd  chord  c(y)  , we  can  write 
down  the  horizontal  separation  x^(y)  of  Tq  and  T in  the  reference 
configuration; 


Xi(y)  = x^(y)  + c(y)  - x^(0)  - c(0) 


We  also  have  the  angle  6 between  T^T^  and  the  x-axis,  defined  by 


tan  0 = Z.J./X.J, 


At  incidence  a , the  line  T„T  rotates  to  T„T  where  T is  the  new  projec- 

0 y 0 

tion  of  the  local  trailing-edge  on  y = 0 . Then  the  total  ground  clearance 
g.j,(y,a)  of  the  local  trailing-edge  is  the  algebraic  sum  of  the  root  ground 
clearance  and  the  vertical  separation  between  Tq  and  T : 


(*T  * ^t/  “ 


h_  + I x„  + z_ 
T V T T 


in  (0  - a) 


■ h.j,  + z^(y)  cos  a - 3^(y)  sin  a 


g.j(y.a) 
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SYMBOLS 


c(y) 

‘^LL<y> 

gjCy.a) 

R . R„ 


chord 

two-dimensional  section  lift  coefficient 
local  section  lift  coefficient 

hei^t  above  ground  of  trailing-edge  of  wing  section  in 
y constant  at  incidence  a 

height  above  ground  of  wing-root  trailing-edge 
height  above  ground  of  quarter-chord  point 
number  of  iterations 

residual  errors  in  symmetrical  and  antisymmetrical  wing  boundary 
conditions 

interference  velocity  (induced  on  wing  by  ground  image 
singularities) 


w 

G 

components 

of 

Hg 

in  (x,y,z)-axes 

w' 

G 

components 

of 

in  (x’,y,z')-axes 

U,  V,  W total -velocity  components 

X,  y,  z axes  in  plane  of  wing  chordline 

x',  y,  z'  axes  in  plane  of  image  chordline 

X,  y,  z free-stream  axes,  centred  on  wing-root  trailing-edge 

Xj^(y)  leading-edge  ordinate 

z wing  thickness 


a 

a* 


angle  of  incidence 
local  twist  angle 

a + 
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Fig. 2 The  local  trailing  - edge  ground  clearance 


Variation  of  lift  coefficient  with  height  above  ground 
at  fixed  incidence  lOVo  RAE  100  section 


Fig. 5a  Effect  of  height  above  ground  on  chordwise 
pressure  distribution:! a)  ti  = 0.  'Airbus’ 
type  wing  (Var.l  ),a=-0  05* 
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Fig. 8 Effect  of  different  trailing-edge  height  configurations 

on  spanwise  centre  of  pressure  locus 
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